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The hydro dynamic equations for a crystals with intersti- 
tials, taking into account the dissipative processes of the vis- 
cosity, heat conduction and the interstitial diffusion are de- 
rived. To achieve that we use the phenomenological approach 
has originally been applied in the theory of the superfluids 
for the derivation of equations of the two-fluid hydrodynam- 
ics. On the basis of obtained equations the problem of the 
propagation of plane waves in the crystal with low interstitial 
concentration has been considered. For the case when effects 
of the viscosity and the heat conduction are absent and the 
diffusion mobility of interstitials is small the absorption coef- 
ficient of longitudinal sound wave has been calculated. 

PACS numbers: 46.05. +b, 62.30.-f-d 



A great number of various in its nature processes occur- 
ing in the bulk and on the surface of the crystalline solids 
can be considered within the scope of the hydrodynamic 
descriptionEl^tiJ. The presence of the slowly varying spa- 
tial and temporal disturbances in the system allows us to 
describe the nonequilibrium behavior in terms of a few 
slow, hydrodynamic variables. One of the most impor- 
tant hydrodynamic processes is the diffusion of intersti- 
tials in the crystalline lattice. Within the scope of the hy- 
drodynamic approximation their dynamics is commonly 
determined by the well known diffusion equation rested 
on the classical Pick's laws. This equation is closed with 
respect to the interstitial concentration and the presence 
of the lattice is only manifested at the calculation of the 
diffusion coefScient for the specific microscopic model. 
Such an approach, justified in the case of the low impu- 
rity concentrations and weak inhomogeneity in the sys- 
tem proves to be unfit at the finite concentration and the 
large gradients. In the case the interaction between the 
impurity system and the lattice is also manifested at the 
macroscopic level through the deformation of the crystal. 
Therefore at the finite interstitial concentration the dy- 
namics of interstitials is to be considered in common with 
motion of the lattice. As a result the hydrodynamic equa- 
tions describing the corresponding process are to contain 
both the interstitial and lattice variables. 

In this paper we shall obtain hydrodynamic equations 
for the crystal containing interstitials. With that end 
in view we use the method haa_|QLiginally been applied 
in the theory of the superfluidsoEj for the derivation of 
the equations of the twp-fluid hydrodynamics and then 
employed for supersolidstil. It gives a possibility to obtain 
equations describing the dynamics of interstitials relative 
to the moving lattice. 

On the basis of the found equations we shall consider 



the propagation of plane waves in the crystal and in 
the simplest approximation calculate the sound veloci- 
ties and absorption coefhcient caused by a diffusion of 
interstitials. 

We assume that the concentration of vacancies is negli- 
gible and the nonequilibrium state of the crystal is char- 
acterized by the mass densities of impurity particles pp 
and lattice atoms p^ and denote by j the current density 
of the medium (lattice plus interstitials). 
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where Vj, and Vp are the velocity fields of the lattice and 
the impurity system, respectively. 

The complete set of the hydrodynamic equations has to 
contain the local conservation laws of mass, momentum 
and equation for the entropy, following from the second 
law of the thermodynamics. We write down the corre- 
sponding balance equations as 
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here p = pL+ pp is the total mass density of the medium, 

V = p^^j is the mass velocity, S is the entropy density of 
the medium, 11^^ is the tensor momentum current density, 
q is the heat current density, J is the diffusion current 
density of interstitials defined as 

J = Pp(vp - v), 

V is the chemical potential of interstitials per unite vol- 
ume, T is absolute temperature and R {R > 0) is the 
dissipative function of the medium; summation over the 
repeated indexes is implied. 

In addition, the system of equations (g) - (Q) has to be 
supplemented by the continuity equation for the intersti- 
tial density and equation of motion for the lattice 
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where c = Pp/ p is the interstitial concentration, f is the 
mass force with which the impurity system acts on the 
lattice; W-Lik is the tensor momentum current density of 
the lattice which we take in the form 



n^ifc = PLVLiVLk - cfik - (y'ik : 



(7) 



where aik is the symmetric, elastic stress tensor of the 
lattice 
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and al describes the effect of the viscosity. Both aik and 
a'^f^ are supposed to be known. Here and further the tilde 
is used to denote the traceless part of a tensor. 

Now let us introduce the vector u defining displace- 
ment of lattice sites and connected with velocity field v^^ 
by relationship 
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The variables p, S, T, c, j, u completely define the 
nonequilibrium state of the-system and satisfy equations 
(0) - (ph . Similarly to Ref .E3 we shall find the remaining 
unknown values 11^^, i?, q, J, f so that the conservation 
energy law 



— +divQ = 0, 
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with E and Q being energy density and energy current 
density, respectively, would follow from Eqs.(y) - (|6|). In 
the Eq.(nQ) the current Q is originally unknown as well. 
To determine the form of the unknown values indicated 
above we shall pass to the new frame moving with veloc- 
ity -Vl , in which the velocity of the lattice of the given ele- 
ment of the medium is equal zero. The energy E and the 
momentum density j are related by the Gallilean trans- 
formation to its values Eq and jo in the frame where the 
lattice rests by relationships 
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Let us write down the differential of _Eo, considered as 
a function of iS, p, c, jo and the infinitesimal strain tensor 



Uij, in the form 



dE^ = TdS + pdp + vdc - 
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wdjo , 



(13) 



here p is the chemical potential of the medium and 
w = Vp — Vj^ is the relative velocity. In the Eq. ( p^ ) 
it has been taken into account that the variation duu of 
sum of diagonal components of the strain tensor Uij is 
determined by the variation of the density dp. 

Differentiating the Eq.(ll) with respect to time and 
using Eqs. (g) - (O), one obtains 
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The last equality after long and tedious, though simple, 
transformations is led to the form 
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where iTik is defined from the equality 

Ilifc = PLVLiVLk + PpVpiVpk +pSik 



TTifc 



and the notation has been introduced 

p = -Eo + TS + pp + w -20 + (JijUij 



(14) 
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In the derivatiopn of Eq. (11^ we have neglected the term 
dijUikdvLk/dxj that is small in framework of the linear 
elasticity theory. In addition, for the terms of the second 
order of infinitesimal in strains the approximation has 
been used 

(vp - Vi) • V{aikUik) ~ Vp • V{aikUik) , 

because the interstitial velocity is appreciably more than 
the lattice one. 

The Eq.(ffl5) defines the momentum current density of 
the medium. The comparison of the Eq.(|l4|) with the 
energy conservation law dTfl) leads to the definition of: 
the energy current Q 
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the dissipative function R 
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and the force f 
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In the framework of the hnear theory the positive defi- 
niteness of R leads to the hnear relationships relating the 
dissipative currents to the thermodynamic forces. Tak- 
ing into account Onsager's reciprocity relations for the 
transport coefficients and time-reversal property of dissi- 
pative effectall, we can write this relationships, at given 
(t'j,, in the form 
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where the fourth-rank tensor 77 is related to the impu- 
rity viscosity, the second-rank tensor k has meaning of 
the pure heat conductivity, /3 is the second-rank tensor 
related to the interstitial diffusion and a is related to the 
cross effect of thermal diffusion of interstitials. 



It is seen from the Eq. ( |15[ ) that p, given by the Eq. ( 16 ) , 
can be interpreted as "pressure" in the medium. The 
relation similar to the Gibbs-Duhem relation for fluid 
system follows from Eq.dlQ). 

pdp = dp + lyd c — SdT — jo^w — Uijdatj . (21) 

From ( pl| ) one obtains for gradients 

pV/i -I- ^Vr = Vp + j^Vc — joiS/wi — Uij'S/aij . 



Introducing this into the Eq.(|19D, we obtain the expres- 
sion for the mass force as 



i = Vp^-{l~c)Vp + cp^V- 



cUijVa-ij . (22) 



This expression shows explicitly that in the limit c — > 
the mass force vanishes provided that p -^ pj^. From this 
condition it follows that it must be 
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The last relationship should be considered as the def- 
inition of the "pressure" pj^ . Let us Jiote that the 
similar expression has been used in Ref.a for a crystal 
with vacancies. In the same limit n^fc coincides with 
the normal expression for a momentum current density 
Ilifc = PLVLiVLk - cTik - (r'iu and in the linear approxi- 
mation in strains the expression dT^) is reduced to the 
standard definition of the energy current density in the 
viscoelastic medium 
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Having substituted the found expression for 11^^ and f in 
(y) and (^ and restricting to the linear terms in lattice 
strains, we obtain the complete set of the hydrodynamic 
equations for a crystal with interstitials 
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(23) 



9j 

— + v^divj + (j • V)Vi + jodivvp 

+ (Vp • V)jo = - Vp + V • CT + V ■ cr' 

-V-TT, 

where R and the dissipative currents are defined by ex- 
pressions (llq) and (00). In addition, the system (p3) has 
to be still supplemented by the conservation energy law 

Now as a simplest example of the application of the ob- 
tained equations we shall consider the propagation of a 
sound in the infinite cristal with the low interstitial con- 
centration. We shall assume that effects of the viscosity 
and the heat conductivity are absent so that the only 
dissipative process is related to the interstitial diffusion. 
In addition, we shall use the isotropic approximation for 
both the elastic stress tensor and the diffusion current, 
taking 



= 2^oUii 
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where p^ is the shear modulus and the constant (3 is pro- 
portianal to the diffusion coefficient. Further we shall 
confine ourselves to the case of the small diffusion mobil- 
ity of interstitials when /3 <C 1. Assuming that a plane 
wave propagates along the a;-direction and having elim- 
inated the current j, one writes the set of the linearized 
equations (E3h in the form 



,2,V 



P 



pc = r/3 
■I 

ps - 
PlU 



dp 
a?" 3 



4 d^u^ 
P-o- 



dx^ 



dx'^ 
V d'^V' 



p dx"^ 



{V = vjpT) , 
{s = SIp) , 



dp' 



92 



(^~^)a^ + 3^°"fe2 

■,d'^u. 



(24) 



2^X 
""' dx^ ' 



dx"^ 



Here the prime denotes small deviations of a variable 
from its equilibrium value which is without prime and 



s — S/p is the entropy per unit mass. It follows from the 
last two equations that in the assumed approximation the 
sound velocity of the transverse waves c^ = ij-o/Pl does 
not depend on the interstitial concentration and coincides 
with its usual value in the isotropic elastic medium. 

The equality (|2l| ) shows that p, c, T, w'^ and aij can 
be taken as the independent variables. Taking into ac- 
count the fact that the scalar function can depend on 
the tensor aij by convolution dijdij only, in the linear 
approxomation, one can consider p, F, s to be functions 
of p, c, and T. In the following we shall take the space 
and time dependence of all variables to be of the form 
exp[iuj[t + X / v)]^ where v is the longitudinal sound veloc- 
ity , and write down the set (EJ) as 
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(25) 



We shall find the sound velocity in the form v = vq + (3vi, 
where /? ^ 0. Having put in (|2^) P ~ one finds zero 
approximation vq as 
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where 0{c) denotes the small terms to be proportianal 
a concentration c. Using the well known properties of a 
functional determinants, one has 

d{p,s) ^ d{p, s) djp, s) ^ /dp\ /ds_\ 

d{p,T) d{p,s)dip,T) KdpJsKdTJp' ^ ' 

Since the adiabatic compression modulus per unit mass 
Kad is defined as 



Kad p\dp 



we have, instead of ( p7| ) 
d{p,s) 



_P_fds_\ 
Kad^dTJp- 



dip,T) 
Substitution of the last equality in the Eq. ( pq ) yields 



2 4/Xo Kad 

3 Pi Pl 



(28) 



since p ^ pi, at c — > 0. The expression (|2^) coincides 
with usual value of the longitudinal sound velocity in the 
isotropic elastic medium. 

The compatibility condition of the set of Eqs.(P5|) up 
to terms of the second order in /3 yields an equation for 
vi . To avoid too combersome expressions we use a num- 
ber of simplifying assumptions. It is known that when 
introducing an impure particle into a perfect crystal its 
volume changes in a macroscopic value. Therefore one 
can suppose that at c ^- dp/dc has to be considerably 
more than all the other thermodynamic derivatives in- 
cluded into the set (p5|). In addition, we shall confine 
ourselves to the temperature region in which one can ig- 
nore a heat expansion and neglect by terms containing 
dp/dT. Then keeping in an expression for vi only terms 
to be proportianal to dp/dc^ we obtain 
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where 
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Allowing for w = uq 
coefficient 



/3fi, one has for the absorption 
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As it seen from Eq.(29) 7; to be proportianal to cj^ and 
the diffusion coefficient. 

In conclusion, we have obtained the complete system 
of the hydrodynamic equations for a crystal with intersti- 
tials. These equations allow to describe the joint dynam- 
ics of the lattice and the impurity system and take into 
account the dissipative efects of the viscosity, the heat 
conductivity and the diffusion of interstitials. On the ba- 
sis of found equations we have considered the problem of 
the propagation of plane waves in the crystal with small 
interstitial concentration and calculated the sound veloc- 
ities and absorption coefficient provided that the viscos- 
ity and the heat conductivity are absent and the diffusion 
mobility of interstitials is small. 



Pershan, Phys. Rev. A6, 



* E-mail: buchb@univer.omsk.su 
^ P.C. Martin, O. Parodi and P.S. 

2401 (1972). 
2 P.D. Fleming and C. Cohen, Phys. Rev. B13, 500 (1976). 
^ B.M. Aizenbud and N.D. Gershon, Physica AIO8, 583 



(1981). 



* T.C. Lubensky, S. Ramaswamy and J. Toner, Phys. Rev. 

B32, 7444 (1985). 
5 A.M. Anile and S. Pennisi, Phys. Rev. B46, 13186 (1992). 
® Z.W. Gortel and L.A. Turski, Phys. Rev. B45, 9389 (1992). 
^ A.M. Anile and O. Muscato, Phys. Rev. B51, 16728 (1995). 
® H.T.C. Stoof, K. Mullen, M. Wallin and S.M. Girvin, Phys. 

Rev. B53, 5670 (1996). 
® M.W. Wu, H.L. Cui, and N.J.M. Horing, Phys. Rev. B54, 

2351 (1996). 
"' E.B. Sonin and W.F. Vinen, J. Phys.:Condens. Matter, 10, 

2191 (1998). 
" I. Tokatly and O. Pankratov, Phys. Rev. B60, 15550 

(1999). 
^^ L.D. Landau and E.M. Lifshitz, Mechanics of Fluids, 2nd. 
ed. Course of Theoretical Physics Vol.6 (Pergamon Press, 
New York, 1987). 
^^ I.M. Khalatnikov, Theory of Superfluids (Nauka, Moscow, 

1971). 
" A.F. Andreev and L.M. Lifshitz, Sov. Phys. JETP 29, 1107 
(1969). 



